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Abstract 

We present a complete analytical resolution of the one dimensional Burgers equation with the elastic 
forcing term —ti^x + /(t), k G M. Two methods existing for the case k = are adapted and generalized 
using variable and function transformations, valid for all values of space an time. The emergence of a 
Fokker-Planck equation in the method allows to connect a fluid model, depicted by the Burgers equation, 
with an Ornstein-Uhlenbeck process. 
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I. INTRODUCTION 



Burgers equation is known to have a lot in common with the Navier-Stokes equation. In par- 
ticular it presents the same kind of advective nonlinearity, and a Reynolds number may be defined 
from the diffusion term In addition, this nonlinear equation is much used as model for statis- 
tical theories of turbulence from which asymptotical behaviours may be determined. But, from an 
analytical point of view, the inhomogeneous form is poor studied, the complete analytic solution 
being closely dependent of the form of the forcing term. For example, the solution of the one 
dimensional Burgers equation with a time-dependent forcing term 

dtu + ud^u - V d^^u = f{t) 
u{x, 0) = (p{x), 

may be obtained by two methods. The first method lies on the Orlowsky-Sobczyk transforma- 
tions (OS) |j2|, where the inhomogeneous Burgers equation ([T]) is transformed into a homogeneous 
Burgers equation. Nevertheless, there exists an other equivalent method to solve analytically this 
problem. By the way of the well-known Hopf-Cole transformation L3J, an inhomogeneous Burgers 
equation may be transformed into a linear equation: the heat equation with a source term, which 
is nothing but a Schrodinger equation with an imaginary time, and a space and time dependent 
potential. Then, several methods have been developed over past decades to treat this kind of equa- 
tions. One of them, the "Time-Space Transformation method" (TST), has been used in order to 
solve the Schrodinger equation with a time dependent mass moving in a time dependent linear 
potential (M. Feng ]^). It is thus shown, ref.lO], the equivalence between the TST method and the 
Orlowsky-Sobczyk method, that is to say, the possibility to solve analytically by two equivalent 
ways the Burgers equation with a forcing term in f(t). The following diagram resumes this equiv- 
alence, where Heat-S designs the heat equation with a source term, BE the Burgers equation, and 
HC the Hopf-Cole transformation. 

OS 

Inhomogeneous BE : f(t) > Homogeneous BE 



HC 



HC 



Heat — S (linear) '^^'^ > Heat 

This yields to present this paper as a continuation of the previous existing methods. The two 
latest methods are adapted in order to solve the inhomogeneous Burgers equation with a forcing 
term of the form —k^x + f{t), where the value represents the string constant of an elastic 



force. Let us note that Wospakrik and Zen have treated this problem but only in the limiting 
case where the diffusion coefficient tends to zero for the asymptotic mode, whereas the methods 
presented here are valid in all cases. The outline of the paper will be thus as follows: the next 
section is devoted to the treatment of an elastic term, firstly by the way of a TST method, and then 
by using a generalized OS method. It is then shown that a Fokker-Planck equation, associated 
to the Ornstein-Uhlenbeck process, arises in the resolution by the TST method. Consequently, 
an "adapted" Hopf-Cole transformation may be obtained for this case, which allows physical 
interpretation in the asymptotic limit. 



n. RESOLUTION FOR AN ELASTIC FORCING TERM 

As underlined in the introduction, the TST method allows to solve a Schrodinger equation for 
some kinds of potentials. So the inhomogeneous Burgers equation has first to be transformed into 
such an equation. Starting from the following one dimensional Burgers equation with a linear 
forcing term 

dtu + udxU — V dxxU = —K^x + f{t) 

(2) 

0) = '•p{x), 

we apply a Hopf-Cole transformation of the form u{x, t) = —2i':^^^dx^{x, t) to obtain a heat 
equation with a source term S: 

dt'^{x,t) = udxx'^ix,t) + S{x,t)^, (3) 

where S{x, t) = — ^x + c{t), c{t) being an arbitrary time-dependent function. This kind 
of equation permits to apply a TST method based on several change of variables. In [5], and 
following [|4|], a TST method has been used in order to solve a Schrodinger equation with a linear 
potential. Here, a quadratic potential appears in Eq. Q, so the method will consist this time to put 

^(x,t) = P(x,t)e'*(^'*), (4) 

with h{x,t) = aix^ + a2{t)x + a3{t) ; ai, a2{t) and a^^t) being constant or time-dependent 
functions to be determined. The transformation ^ introduced in Eq. Q gives 

dtP = V dxxP + 2iy dxh dxP + (z/ dxxh + p^d^hf + S- dthj P. (5) 
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Then, in order to cancel the factor of P, we put 



ud^^h + v{d^hf + S-dth = 



(6) 



which gives a polynomial of second degree in x. This polynomial becomes zero since all its 
coefficients are. It comes respectively 

,^2 



Auaf + — = 0, 
/ 

4z/aia2 (22 = 0, 

2z/ai + ua^ + c — = 0. 



(7a) 
(7b) 
(7c) 



Since Eqs. O are satisfied, Eq. ^ is simplified to 



dtP = V d^^P + 2z/ d^h d^P. 



(8) 



We now apply to Eq. ^ the following change of variables 

y = r{t)x + 
t' = t. 

This induces a transformation of Eq. (EJ) into : 



(9) 



dt,P =Ur^dyyP + 



-r jr + 4i^ai)(?/ — g) + 2ura2 — g 



We have now to cancel the term in dyP, so we put 



(10) 



r — Auair = 0, 
2vra2 — g = 0. 



(11a) 
(lib) 



Notice that the relation (TTal) gives 



where i= v — 1, with the result that the solution of Eq. (II lal) will be 



(12) 



(13) 



Eqs. (fTTT) being satisfied, we obtain 



i9t'P = i^r^dyyP ; 



(14) 



and finally the transformation 

T{t') = [ r\s)ds , (15) 



yields to the expected heat equation: 

drP{y,r) = udyyP{y,r). (16) 

We show now that the Orlowsky-Sobczyk method is a particular case of the method employed 
here for an elastic term: the Generalized Orlowsky-Sobczyk method (GOS). 
Let us consider again Eq. Q, and let us introduce a new velocity v = v{x, t) such as 

u = vr{t) +ax + ijj{t) , (17) 

where r(t), a, ipit) are time dependent functions or constant determined later. The transformation 
dm introduced in Eq. © yields to : 

V (r + ar) + x^K^ + a;^j + ^■j/' + a^/' — +rdtV + r'^vdxV + arxdxV + ripdxV — urdxxV = 0. (18) 

In order to delete the terms in v and x, and those only depending on time, we put 

r + ar = (19a) 
+ = (19b) 
ip + aip- f = (19c) 

Since the system (fT^ is verified, then Eq. ( fTSl is simplified into 

rdtV + r'^vdxV + arxdxV + ripdxV — urdxxV = 0. (20) 
Then, the same time and space change of variables as Eq. ® applied to Eq. (I^DI) leads to 

pdt'V + (rg + r'^ip) dyV + (r + ar){y — q)dyV + r'^vdyV — ur^dyyV = 0. (21) 
After what, putting 



rq + r^?/^ = (22) 



we obtain 



If we put now t' as 



\dt'V + vdyV = vdyyV . (23) 



t' 

T{t') = I r\s)ds , (24) 







it comes a homogeneous Burgers equation governing the new velocity v : 

drV + VdyV = V dyyV . (25) 

From this, the HC transformation v = —2v^dyP yields again to the expected heat equation 

drP{y,r) = udyyP{y,r). (26) 

Hence, both methods GOS and TST may be connected thanks to a commutative diagram like the 
one of the introduction, with a force —k^x + f{t). 



III. DERIVATION OF AN ORNSTEIN-UHLENBECK PROCESS 



Let x{t) be a stochastic variable satisfying the following Langevin equation and describing an 
Omstein-Uhlenbeck process 



— — = —Kx + y2uhit)\ 



(27) 



where h{t) stands for a Gaussian white noise verifying the standard conditions 

{h{t)) = and {b{t)b{t')) = 5{t - t'). 



(28) 



Then, using a Kramers-Moyal expansion, a Fokker-Planck equation may be obtained for the tran- 
sition probability P(x, t) 



dtP{x,t) = {xP{x,t)) + z/9^^P(x,t). 



(29) 



This equation is usually solved by Fourier transform, and the solution P = P(x,x',t) for the 
initial condition P{x, t\x', 0) = 5{x — x') reads 

K,{x — e~'^*x'Y 



P 



K 



27Tu(l- e 



exp 



2u{l 



) 



(30) 



It is shown in appendix that this solution may also be found by the TST method. 
The interesting point lies in a connexion existing between the Ornstein-Uhlenbeck process (Eq. 
( l29b ) and the Burgers equation Q with f{t) = 0. In order to see this fact, we apply the transfor- 
mation 

P{x,t) = '^{x,t)e-^, (31) 



6 



to the Fokker- Planck equation d29t , which leads to the heat equation 



2 2 
r\j r\j Jb 



dt<^ = vd,,^^ + - ^. (32) 

So, the Hopf-Cole transfomiation 

u{x, t) = -2u^^d,^{x, t), (33) 
transforms Eq. (B^ into the inhomogeneous Burgers equation 

dfU + udxU = I'dxxU — K^x. (34) 

This interesting result implies two remarks. Firstly, this connection gives rise to a physical mean- 
ing of the TST method. Indeed, the function P introduced in the transformation Q is no more 
an unspecified variable, but takes the sense of a transition probability for the variable x{t). Then, 
considering both Eqs. (OTT) and (B^ . we obtain a relation between the velocity u and the transition 
probability P: 

u{x, t) = -2u^^—^dxP{x, t) - Kx, (35) 

which is composed of a Hopf-Cole part and of a linear part. So, this relation may be considered as 
a Hopf-Cole transformation adapted to the Ornstein-Uhlenbeck process. Moreover, the asymptotic 
limit of P(x, x', t) is given by (l30b : 

lim P(x, x', t) = J - — exp — — , (36) 
t^oo y 27rz/ \ Iv J 

and thus, from the relation (B51) . we can see that the asymptotic limit of the velocity will read 

lim u{x,t) = Kx, (37) 

which is a stationary solution. The initial condition P{x, t\x', 0) = S{x — x') expressing the fact 
that a particle cannot be at several positions at the same time, it may be considered as the more 
acceptable condition for P. Then, the asymptotic solution (l37t have a real physical sense. We can 
conclude on the fact that an elastic forcing term applied to the system gives rise to a stationary 
transition probability in the asymptotic mode. Consequently, the effects of the oscillations will 
decrease, up to disappear in the long time limit, and stabilize the system with a velocity propor- 
tional to the displacement. The evanescence of the effect of the force is due to the initial condition 
sensitivity of the Burgers equation. We can see thereby on the system, a phenomenon closely 
connected to the turbulence effect: the lost of memory in the long-time limit. 
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IV. CONCLUSION 



We have presented the complete analytical solution of the Burgers equation with an elastic 
forcing term. The methods presented here have been used before but only in the case of a time- 
dependent forcing term. As a perspective, we can say that the generalisation of the methods to any 
order of power of x seems actually be a difficult task. Indeed, a transformation of the form y 
r{t)x + q{t), has been introduced in order to delete terms proportional to x. So this transformation 
seems without effect when higher powers of x appear. Moreover, the more the degree will be 
high, the more the resolution will be difficult, due to the increasing number of variables to be 
introduced. The second main result of the paper lies in the existence of links between a fluid 
model (Burgers) and the statistical physics (Ornstein-Uhlenbeck). By a set of transformations, 
we have connected the Burgers equation for the velocity u = ds/dt to a Fokker- Planck equation 
for the transition probability of the variable x. From the Burgers equation (l34b . the transformation 
allows to get directly the Fokker-Planck equation ( E^ as a specific Hopf-Cole transformation. 
It appears that the linear force, describing the Omstein-uhlenbeck process, stabilize the system in 
the asymptotic mode with a velocity proportional to the force applied initially, since we consider 
the initial condition P{x, t\x', 0) = S{x — x') as the more acceptable condition. This result shows 
a characteristic property of turbulence, i.e the unpredictability in the long time limit of a velocity 
field governed by the Burgers equation. An application of the methods presented here will be 
described in a forthcoming paper with the case of an electric field in a plasma. 
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APPENDIX: SOLUTION OF THE ORNSTEIN-UHLENBECK PROCESS 



We show that we can recover the solution (BUb by the way of our TST method. 
Rewriting Eq. (1^ . 

dtP = I'dxxP + KxdxP + K-P, 



we apply the change of variable 



y = r{t)x, 
t' = t. 



This yields to 



dt,P = ur^dyyP + U - - ydyP + nP. 



To cancel the term in dyP we put obviously 



ft; - 1 = ^ rit') = e"^' 
r 



This leads to 



Then, putting 



followed by the transformation 



we obtain the heat equation 



dt,P = ur^dyyP + KP. 



Piy,t') = e{y,t')e^'\ 



Tit') = / r'is)ds, 



drO = vdyyQ. 



Notice that the condition P{y, y', 0) = 5{y — y') implies that 6(?/, y\ 0) = 5{y — y') 
mental solution of (IA.8I) is thus 



e(y,r) 



exp 



TIVT 



r\2 



jy - y' ) 



(A.l) 
(A.2) 
(A.3) 

(A.4) 

(A.5) 

(A.6) 

(A.7) 

(A.8) 
The funda- 

(A.9) 



after what, putting y and r in place of their expression, it is to say 



y = xe , 



r = k (e^''*' - 1) 



we obtain 



P 



K 



-2Kf 



exp 



— e ^^x'Y 
2z/(l - 6-2"*) 



(A. 10) 



(A. 11) 
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which is the same result as Eq. (l30l) . 
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